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1. Introduction 



1.1. Polytopes and Lattices 

Our main goal is the classification of smooth lattice polytopes with few lattice points, 
so we will start by introdncing polytopes and lattice polytopes. 

Definition. A set P C M'^ is called a polytope if it is the convex hull of a finite set of 
points V = {f 1, . . . , Vn} C M"*. If V is minimal up to inclusion, then the points Vi, . . . ,Vn 
are the vertices of P. 



Remark. In the following we will always assume V to be the vertex set of the polytope. 

For the definition of smooth lattice polytopes, and to make precise what we mean by 
polytopes with few lattice points, we first need to define lattices in M.'^. 

Definition. For some finite set B = {bi, . . . , bn} C M*^, we define the lattice A.{B): 



If the vectors in B are hnearly independent then B is called a basis of A(S). 

Remark. We will now restrict to the lattice A — Z^C M.^, so any set B — {bi, . . . , b^} C 
Z'^ is a basis for Z'^ if det(S) = ±1. 

Now, for a lattice A = A{B) C W^, the dual lattice is 

A* = {7 e (R'^)*| 7(x) e Z Vx e A} 

and one can easily see that for A{B) — Z'^ we have A* = A{{B~^)'^) = Z'^. So for some 
lattice basis B define the dual lattice basis B* — {B~^Y . 

We can now define smooth lattice polytopes: 

Definition. Let P be a polytope. Then P is called a lattice polytope if all vertices 
have integer coordinates, i.e. V C Z"^. 

When we talk about the lattice points L{P) of a polytope P C M*^ we mean the elements 
of Z"^ that lie in the polytope, namely L[P) — P n Z^. So we want to limit the number 
|L(P)| for our polytopes. 

For a lattice polytope P to be smooth it has to be simple. So for v denote by nb(f ) 
the set of adjacent vertices, i.e. for every w G nb(?;) there is an edge containing v and 
w. Then a polytope P C M*^ is simple if | nb(t')| = d for all vertices f G V. 
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Definition. A simple lattice polytope P is called smooth if, at each vertex, the lattice 
minimal edge-directions form a lattice basis, or more formally 



1.2. Cones and Fans 

In our main algorithm we want to approach the polytopes from the dual and therefore 
we now introduce cones and fans and define some important properties. 

Definition. For S finite, the polyhedral cone C — cone(5') is 



A cone C C M*^ is called pointed if there is a G (R'')* such that 

Cn{a; e E'^ I a{x) < 0} = {0} 

If C = cone(5') is pointed and S is minimal up to inclusion, then the vectors r & S are 
called the rays of C. 

Remark. Some restrictions on the cones we will consider later on. 

• All our cones are pointed. 

• All our cones are rational, which means that there is C Q*^ finite, such that 
C = cone(5). 

• And we assume all rays r e 5" to be primitive vectors, i.e. gcd(r) = 1 for all r & S. 

The dimension of a cone C = cone(S') is the dimension of its linear hull dim(C) = 
dim(span(S')). A A;-dimensional cone C is called simplicial if there is S* = {ri, . . . ,rfc} 
such that C = cone(5'). One may already see the corresondence to simple polytopes. 
We will continue in this direction by defining smooth cones. Later we will see that both 
definitions of smoothness coincide if we choose the right construction to define a fan 
from a polytope. 

Definition. Let C C K'' be a /c-dimensional simplicial cone, then C is smooth if there 
are Zi, . . . , z^. E Z'^ such that 

• C = cone{{zi, . . . , Zk}) 

• det{zi, . . . , Zd) ^ ±1 




conei 




ttrT e W\ Qr G R>o E S 



4 



Definition. A fan F is a finite set of cones witli tlie following properties. 

(1) Ci,C2GF^CinC2GF 

(2) C eF, 0, C <C (i.e. C is a non-empty face of C) ^ C e F 

For a fan F, denote by F^''^ the set of cones of dimension k. A fan F in R'' is called 
complete if it covers the whole M*^, which means UceF ^ ~ 

The properties of cones defined above can of course be extended to fans, a fan is sim- 
plicial or smooth if all its cones are. 

1.3. Duality 

In the third section of the introduction we will now state a duality theorem of polyhedral 
theory and establish a connection between polytopes and fans. This will turn out useful 
later in the classification. 

Definition. For a matrix A G M"^^"^ and a vector b G the polyhedron P{A, b) C 
is the intersection of finitely many halfspaces 

P{A, b)= Pi {x G W^l {au x) < bi} = {xe R'^\ Ax < b} 

CL'l row of A. 

For the duality theorem we need the Minkowski sum P + Q oi two subsets P,Q G M'^: 

P + Q = {x + yeR'^\xeP, yeQ} 

Theorem. A set P C M'^ is a polyhedron if and only if it is the Minkowski sum of a 
polytope and a cone. 

Proof. The proof of this duality can be found in most textbooks containing polyhedral 
theory, e.g. Lectures on Polytopes [8]. 

Remark. In particular we have the following results for polytopes: 

• A polytope is a bounded polyhedron. 

• For P = P{A, b) bounded, there is ^ C M"^ finite, such that P = conv(l^). 

• And for P = conv(\/) there are A G W^""^ and 6 G M™ such that P = P{A, b). 

This is very useful as it allows us to switch between both representations as we need it. 
Now let us define the normal fan of a polytope, the main construction we will use for 
the classification of the polytopes. 
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Definition. Let P — P{A, h) C R*^ be a polytope, the normal cone of a face F c P is 

NCp(F) = cone({ai row of A \ {ai,x) = hi Vx e F}) 
The normal fan of P is defined by the normal cones of all faces of the polytope. 

NFp = {NCp(F) \FcP face} 

To sec that this does indeed produce a fan, observe that for every face F of the polytope 
the faces containing F are defined by subsets of the rays defining F. Thus, the faces 
of P containing F define the faces of NCp(F). Additionally this construction always 
produces a complete fan. 

This yields a construction to build a fan from a given polytope. We can 
also use this in the other direction. For a given matrix A containing 
the rays of the fan we can define a polytope P{A, b) by choosing a 
right-hand side vector h. We just need to make sure that the normal 
fan of the polytope is again the fan we started with, since there can be 
different polytopes with the same set of facet normals. 
However, we only want to classify smooth polytopes, so the obvious question is whether 
there is a property of the normal fan that coincides with the smoothness of the polytope. 
And the obvious answer is also correct: 

Proposition. Let P be lattice d-polytope and NFp its normal fan. Then P is smooth if 
and only if NFp is smooth. 

Proof. First we will show that P is simple if and only if NFp is simplicial. 

• Let NFp be a simplicial fan, then the normal cone of a vertex v of P is defined by 
exactly d rays. Each (d — l) subset of the facets defined by these rays can define an 
edge and since there are d such subsets, there can be at most d edges containing v. 
Because P is a rf-dimensional polytope, there must be at least d edges containing 

V and thus P is simple. 

• Let P be a simple polytope and f be a vertex of P. Now, every facet containing 

V is itself a (c? — l)-dimensional polytope and thus there must be at least {d — 1) 
edges containing v in each facet. Since there are only d {d — l)-subsets of the d 
edges at v there can be at most d different facets and thus NFp is simplicial. 

We choose again one vertex v in P. To prove the smoothness we use the matrizes 
E, A defined by the edge-directions ei, . . . , at v and the normal cone of the vertex 
Np{v) = cone({ai, . . . , a^}). 

For given smooth edge-directions E choose A = — (P"^)^, obviously A defines the 
normal cone. By Cramer's rule and since det(P) = ±1 we know that A is integral and 
det{A) = =f1) thus the cone is smooth. And for given smooth normal fan with rays in 
A choose E = By the same argument E e Z'^^'^, det(P) = ^1 and thus E 

contains the lattice minimal edge-directions. □ 
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2. Finiteness 



2.1. Finiteness Theorem 

Now we can state the main theorem [3]. 

Theorem. For integers d and N , there are only finitely many smooth d-polytopes with 
at most N lattice points. 

Since the following proof inspired the technique for the classification algorithm following 
subsequently, we will repeat its details. 

Proof. We will proof this in three steps. First we want to limit the number of combina- 
torial types of smooth fans, then the number of smooth fans with each of those types 
and finally the number of smooth polytopes having one of those as normal fan, all using 
the limit on the number of lattice points. 
In the following we will again assume all rays to be primitive. 

• Assume we have a combinatorial type of a complete simplicial fan F with n rays and 
at most N full-dimensional cones. A full-dimensional cone of this fan can be defined by 
d-subsets C C {1, . . . ,n}, \C\ = d. The fan can then be defined by the set of all such 
full-dimensional cones 

^(d) ^ Note that we use a different notation here since 

we consider only the combinatorial type of the fan and do not care about the rays. We 
know that the fan has at most N full-dimensional cones, so IF^"^)] < A^. Additionally 
every ray needs to appear in at least d full- dimensional cones: 

d, = \{C e F'^'^^ \ i e C}\ > d ViG {!,..., n} 

By counting all appearances of all rays in F^'^^ in two different ways we have: 

n 

d-n<^di<d- N 

i=l 

And thus n < N, which means that we have at most rays in the simplicial fan F. 
Hence our fan is fully defined by choosing 

Since A^ and d are fixed there are only finitely many such subsets, and thus only finitely 
many possible combinatorial types of simplicial fans. 

• Given the combinatorial type we now want to define the fan F up to isomorphism: 
Let C = cone(ni,n2, . . . ,nd-i) G -F be a codimension 1 cone which is contained in two 
full-dimensional cones Ci = cone(C, ri) and C2 = cone(C, r2). Since all the cones are 
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smooth ri and r2 have lattice distance one from C, but he on different sides of C. So 
ri + r2 hes in the span of C and thus 

d-l 

ri + r2 = ^ aciUi 

i=l 

The coefficients ac\i are the edge-paremeters of a fan F. We know ac,i G ^ because C 
is smooth and thus {rii, . . . , Ud-i} is a Hilbert basis of C. 

We can also use the edge-parameters ac,i ior 1 < i < d — 1, Ce F^^~^^ together with the 
combinatorial type to reconstruct fan F using the following technique: Every smooth 
fan can, while preserving the lattice, be transformed to contain the cone defined by all 
d unit vectors. This means wc can assume the fan to contain this cone. Then there 
are codimension 1 cones for every (d — l)-subset of the unit vectors, so, given the edge- 
parameters for each one of those cones, one can calculate the ray on the other side of 
this codimension 1 cone with the formula above. Iterate until all rays are calculated. 
Since the edge-parameters are in Z there may be infinitely many such smooth fans, but 
in the next step we will show a way to bound those parameters. 

• In the last step we have to bound the number of smooth polytopes P with at most 
N lattice points having a given normal fan F and the number of smooth fans the above 
construction yields. 

We therefore consider the edge e of the polytope dual to a codimension 1 cone C E F. 
Let /(e) be the lattice length of the edge e and Xi,X2 be the vertices defining this edge. 
Denote by Ui,...,U(i the lattice-minimal edge-directions at Xi which are dual to the 
inner normals —rii, . . . , —Ud-i, —ri around xi. So we have X2 = xi + l{e)ud. Using the 
edge-directions at xi and the edge-parameters we can calculate the edge-directions at 

X2. 

Lemma. For Xi,X2 as above, the edge-directions at X2 are 

u'i^Ui + aciUd for i e 1} 
u'd = -Ud 

Proof. The inner normals around X2 are —rii, . . . , —nd-i, —r2 and 

d-l 

-r2 = ri-^ aciUi 
1=1 

so we have 

{-rii, u'j) = {-rii, Uj + ac,jUd) = (-rij, Uj) - ac,j{ni, Ud) for i,j <d 

= {-rii, Uj) = 5ij for i,j < d 

{-rii, u'd) = {-rii, -Ud) = ioT i <d 
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d-l 

(-r2, Uj) = {ri,Uj + ac,jUd) - ^ ac,i{ni, Uj + acjUd) for j < d 

i=l 

= acj{ri, Ud) - ac,j{nj, uj) = for j < d 

d-l 

(-r2, u'^) = (ri, -Ud) - ^ ac,i{^i, -Ud) = {-ri,Ud) = 1 

1=1 

which proves that the matrix containing these edge-directions is inverse to the given 
inner normals, finishing the proof of the lemma. □ 

Since in our smooth lattice polytope every edge must have length at least 1 the lattice 
points defined by the vertex X2 plus one of the edge-directions must of course lie in the 
polytope and thus 

{ri,X2 + ac,iUd) > (ri,xi) 
(ri, xi + lie)ud + ac,iUd) > {n, Xi) 
(/(e) + ac,i){ri,Ud) > 

ac,i > -lie) > -N 

This gives a lower bound for the edge-parameters. For an upper bound let us first 
introduce the thickened edge. 



Xi X2 = xi+ l{e)u3 




Figure 1: A thickened edge 

Definition. Using the notation from above we define the thickened edge in a smooth 
polytope as the convex hull of all those points we just discussed. 

c(e) = conv({xi, xi + ui, . . . , xi + Ud-i, 

X2, X2 + U1 + ac,iUd, ...,X2 + Ud-i + acd-iMd}) 

We obviously have c(e) C P and as shown in Figure [H c(e) contains d — 1 lines from 
xi + Ui to X2 + Ui + ac,iUd each containing /(e) -|- 1 -|- ac,i lattice points and the line from 
xi to X2 containing /(e) + 1 lattice points, so 

d-l 

|c(e) n Z'^l = d{l{e) + + 
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Since the edge length must be positive and we want the number of lattice points of P 
and thus of c(e) to be less than N, this gives us an upper bound for the sum of all 
edge-parameters. Together with the lower bound this leaves us with only finitely many 
choices for the edge-parameters. This means for an upper bound on the number of 
lattice points in the polytope there are only finitely many possible smooth normal fans. 
Now for a given normal fan we can define the polytope up to isomorphism by fixing all 
edge lengths. Just assume one vertex at the origin and the edge-directions to be the 
unit vectors, then the edge lengths define the next vertices and the normal fan defines 
the next edge-directions. Iterate until all vertices are fixed. And we want of course the 
edge length /(e) to be positive and less than A^, so for every edge there are only finitely 
many possible edge lengths. And thus only finitely many polytopes for this fan. 

So altogether we have finitely many combinatorial types of fans, finitely many choices 
for the parameters and finitely many choices for the edge lengths and thus only finitely 
many smooth lattice d-polytopes with at most lattice points. This completes the 
proof of the finiteness theorem. □ 

2.2. General Algorithm 

The proof gives rise to the following approach for the classification of such polytopes. 

2.2.1. Generating Normal Fans 

The first step is to find all combinatorial types of complete simplicial fans in dimension 
d with at most full-dimensional cones. One can then proceed by finding parameters 
ac,i which produce smooth fans, for all codimension 1 cones. 

In dimension three we can use that the combinatorial types of simplicial fans are equiv- 
alent to the combinatorial types of triangulations of 5*^ with at most A^ 2-simplices, 
which have been classified by R. Bowen and S. Fisk in Generation of triangulations of 
the sphere [2]. And Tadao Oda has done the classification of smooth two-dimensional 
fans which are minimal up to equivariant blowing-ups and the same in dimension three 
for at most eight rays in Convex Bodies and Algebraic Geometry [7]. We will explain 
and use those results later in chapter 3. 

2.2.2. Generating Polytopes 

The next step now is to find the polytopes. In the proof we used the edge lengths 
together with the normal fan to define the polytope. Here we will use the right-hand 
side of its inequality description P{A, b) = {x \ Ax < b} instead of the edge lengths, but 
after expressing the edge length in terms of the 6-vector we can use the edge lengths to 
find all those 6-vectors defining smooth lattice polytopes with the correct normal fan. 
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Right-Hand Side 



For a given fan F we now need to find all possible 6-vectors which yield smooth lattice 
polytopes with at most lattice points and F as normal fan. Therefore define the set 
of all such vectors 

RHS(F, N) = {beZ"'\ NFp(^,5) = F, \P{A, b)nZ'^\<N} 

where the rows of A are the rays of F. In the following we want to define a polytope 
containing this set and enumerate all its lattice points. 

Calculating the edge lengths 

Lemma. Let C — cone(ni,n2, . . . ,nd-i) be a codimension 1 cone and ri,r2 G 'Z^ such 
that Ci = cone(C, ri), C2 = cone(C, G F are the full-dimensional cones containing 
C. Then the length /(e) of the edge e in the polytope P{A, b) dual to C can be expressed 
in terms of the right-hand side b as follows: 

d-l 

/(e) = br^ + &r2 - X] acibm 

i=l 

where ac,i <i < d — 1) such that ri + r2 = Y^fzl o-c,in'i- 

Proof. The coefficients a* are the same as in the proof of the theorem. Now let Xi and X2 
be the vertices dual to Ci and C2 respectively. Prom the proof of the finiteness theorem 
we know (ua, — ri) = 1 and thus /(e) = {x2 — Xi, — ri). Since ri + r2 — YlfZi (^0,1^^ this 
yields 

/(e) = (xi,ri) - (x2,ri) 

d-l 

= &ri - {X2,^ aciUi - r2) 
1=1 

d-l 

= &ri - X ac,i{x2, rii) + {x2, r2) 

i=l 
d-l 

— Ki ~ X^ 0'C,i^ni + K2 
1=1 

□ 

Approximating the Right-Hand Sides 

Since we can express /(e) linearly in terms of the 6- vector we can now use bounds on the 
edge length to approximate RHS(F, N). For a codimension 1 cone C, denote by ec the 
edge in the polytope dual to C. 
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Lemma. For a given normal fan F in dimension d and a maximal number of lattice points 
A^, every polytope defined by F and a vector from RHS(F, A^) satisfies the following 
inequalities 

Uj 

l{ec) > 1 VC e F('^-i) 

6i = yi<i<d 

{Kec)-l)<N-\F^'^\ 

Using l{e) as defined in tfie previous lemma, these inequalities define a polytope 

B{F,N) D RHS(F,A^) 

Proof. We first show that RHS(F, N) satisfies the inequahties: 

(i) We can now use an idea from the proof of the finiteness theorem. The convex hull 
of the thickened edge c(e) contains d{l{e) + + ^c,i lattice points for an edge 
e with the parameters ac,i, so we have the following bound for l{e). 



(ii) The fan completely defines the combinatorial type of the polytope, so every codi- 
mension 1 cone must define a non-degenerate (i.e. /(e) > 0) edge in the polytope 
and since all vertices have integer coordinates the edge length must be at least 
one. This gives the second set of inequalities. 

(iii) We now assume F to be lattice-transformed such that the first d rays are the 
unit vectors, this is possible because F is smooth. Since we do not care about 
translations, we can fix the first d right-hand sides to be zero. 

(iv) The total number of lattice points on all edges of the polytope is the number of 
relative interior lattice points of each edge (/(e) — 1) plus the number of vertices, 
so this sum must be less than the given bound N on the number of lattice points. 

So all inequalities hold for any b e RHS(F, A^), they are all linear in terms of b and thus 
B{F,N) is a polyhedron. 

Now to prove that this polyhedron is bounded we use a similar approach as we used to 
construct the fan from the edge-parameters. For every codimension 1 cone we have the 
following inequalities 

d-l 

1 < /(e) = bri +br2-Y ac,ibni < N 

1=1 
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Starting with the full-dimensional cone defined by the unit- vectors, where all bi [1 < 
i < d) are zero, we can always find codimension 1 cones where only one of the b^, in the 
inequalities above is not yet bounded and thus bound it. We can iterate this until all 6* 
are bounded and thus B(F, A^) is a polytope. □ 

Enumeration 

Now, given an inequality description of B{F,N) we can calculate all lattice vectors 
b G B(F, A^)nZ™. Together with the matrix A, each one defines a smooth lattice polytope 
P{A, b) with possibly too many lattice points. Thus, after removing all polytopes with 
more than N lattice points, we obtain a list of smooth lattice polytopes with normal fan 
F and at most lattice points. The last step is to remove all polytopes that are lattice- 
isomorphic to another one. All these computations can be done with polymake [6j using 
interfaces to 4ti2 [1], Latte macchiato and normaliz2 [3]. 
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3. Classification 



After a short definition of blowing-ups, the sections 3.2 and 3.3 contain a detailed de- 
scription of how to classify smooth two- and three-dimensional lattice polytopes with at 
most 12 lattice points based on the technique described above. 

3.1. Blowing-ups 

We want to use the results of Tadao Oda who classified smooth fans in dimension two 
and three which are minimal up to equivariant blowing-ups. This means that all normal 
fans we need can be generated by successive equivariant blowing-ups of the fans given in 
Convex Bodies and Algebraic Geometry by Tadao Oda [7]. But first we need to clarify 
what such a blowing-up is. In short, it means subdividing a cone and all cones containing 
it. 

Definition. Let C be a cone in the fan F and r & C. The blowing-up 6(C, r) of a cone 
C at r is 

r) = {cone({r, C'}) \ C E F,C' < C,r i C'} ii r e C 

6(C,r) = {C} ifr^C 

where C < C means that C is a face of C . 
Now the blowing-up r) of F at r is 

h{F,r)= \Jh{C,r) 

C&F 

One can easily see that this is indeed again a fan 

A blowing-up is equivariant if r = Xlsgs ^ some cone C = cone(S') G F, where S is 
minimal up to inclusion and all rays are primitive. We denote by b{F, C) the blowing-up 
in cone C. 

For a smooth fan, an equivariant blowing-up is again smooth. Since we want to classify 
smooth polytopes and thus smooth normal fans, we only have to consider equivariant 
blowing-ups. 

Remark. Since we have a connection between the fan and the polytope one could ask 
what such a blowing-up of the normal fan corresponds to on the polytope side. 
We first choose a cone in the fan, which is equivalent to choosing a face in the polytope, 
and then add the sum of all rays of the cone as new ray. This corresponds to adding a 
new facet by cutting-off that face of the polytope. 
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3.2. Smooth Polytopes in Dimension Two 



3.2.1. Classifying the Normal Fans 

As already mentioned we do not start with combinatorial types of fans, but with minimal 
smooth fans, and in dimension two there are only two such fans as given by the following 
theorem. 



Minimal Fans 

Theorem. (Tadao Oda jT]) Every smooth fan in dimension two can, up to isomor- 
phism, be created by finitely many successive equivariant blowing-ups of one of the fol- 
lowing fans. 

(i) The fan Fp representing the complex projective space P2(C). 

F« = {ao, ai, as} = {(1, 0), (0, 1), (-1, -1)} 

F^''-^ = {cone({ao, ai}), cone({ai, 02}), cone({a2, ao})} 

(a) The fans F^ representing the Hirzebruch surfaces for a E'L. 

= {ao, ai, a2, a,] = {(1, 0), (0, 1), (-1, -a), (0, -1)} 
^i'^^ = {cone({ao, aj), cone({ai, as}), cone({a2, ag}), cone({a3, ao})} 




(-1,-3) 



Tree of Blowing-ups 

Instead of generating combinatorial types or finding parameters we have a different task 
now, namely equivariant blowing-ups. In dimension two this means choosing one full- 
dimensional cone C = cone(ri,r2), inserting the ray r = vi + r2 and replacing the cone 
C with Ci = {ri,r} and C2 = {r,r2}. 

So beginning with a fan F we can generate all normal fans with the following recursion: 
1. For every full-dimensional cone C G F^'^^: 
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a) Create the blowing-up b{F, C), 

b) If b{F, C) has less than 12 full-dimensional cones: Start recursion. 

Which is a depth-first search in a rather large rooted tree, having 1 + (3! + 4! + ■ ■ ■ + 
ll!)/2! = 21977356 nodes for up to 12 full-dimensional cones for the fan Fp, but lots 
or even most of those fans are isomorphic. So we need a way to reduce the number of 
nodes to consider: 

We first assign an order to all full- dimensional cones. Now, since blowing up Cj) in 
Cj and Cj) in Cj will give the same fan, we can ignore all cones Cj with index j < i 
in the subtree of the blowing-up 6(F, Cj). (When blowing up in cone Cj in a fan with k 
full-dimensional cones Ci, . . . , C^ we will index the two new cones with i and A; -|- 1.) 
This reduces the amount of blowing-ups to 58785 for Fp and 35072 for Fa- 

Pcirameters 

The fan Fa and all further blowing-ups of it represent classes of infintely many fans with 
a G Z. But, from the proof of our main theorem, we know that for some upper bound 
N on the number of lattice points, there can be only finitely many of those. The proof 
even gives us a direct way to bound the parameter a. 

Lemma. For the fan Fa and a maximal number of lattice points N for the polytope, we 
can restrict atoQ<a<N and a^l. 

Proof. For a pair of two-dimensional cones Ci = {ri,r} and C2 = {'"2,'"} sharing a 
common ray r we can write ri + r2 in terms of r, i.e. ri + r2 = a^r. Additionally we 
know ar > —l{e), where /(e) is the length of the edge e dual to the ray r. Obviously 
/(e) < N, so we have > —N. 

m For the cones Ci = {(1, 0), (0, 1)}, C2 = {(-1, -a), (0, 1)} this yields 

n + r2 = (1, 0) + (-1, -a) = (0, -a) = -a ■ (0, 1) 

and thus —a — a(o,i),i > —N. 
• And for Ci = {(1, 0), (0, -1)}, C2 = 
n + r2 = (1, 0) + ( 
and a — a(o,-i),i > —N. 

Because of the symmetry of the fan we can restrict to positive a and skip a = 1 since 
Fa then equals Fp after an equivariant blowing-up in the cone C = cone({ao, 02})- □ 



= {(—1, —a), (0, —1)} we have 
-1,-a) = (0, -a) = a - (0,-1) 
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Pcirameters in the Blowing-up 

We could find bounds on the parameter a of the Hirzebruch surface using an edge of 
the corresponding polytope, but when blowing up the fan, which means cutting off a 
vertex of the polytope, this edge may have different parameters and the bounds may 
be different too. We could of course do the same calculation for all blowing-ups, but a 
faster approach is to find a way to use the bounds we found for some fan F on all its 
blowing-ups. 

Assume we have a pair of cones Ci = {ri,r}, C2 = {r2,r} G F around a ray r with 
ri + r2 = ttrV, this means our bound is > —N. Now, in the blowing-up 6(F, C2), we 
have a new cone C2 = {r2, r} with ri^ — r + r2- This yields 

ri + r2 = ttrT — r2 + r + r2 = (or- + 1)?^ 

So we now have ar + 1 > —N and thus > — — 1. This means we have to loosen the 
bounds on a parameter by one if it appears as edge-parameter at one of the rays of the 
current cone. 

Thus, after creating a new fan for every valid value of a for each class of fans, we now 
have all possible smooth normal fans for the smooth polytopes. 

3.2.2. Creating the Polytopes 

Now we can continue as in the general case, for every fan F we define the polytope 
B(F, N) approximating RHS(F, N) and enumerate all its lattice points b. For some fan 
with rays A and a right-hand side b check if \P(A,b) n < N. And finally remove 
polytopes that are lattice-isomorphic to another one. 

3.2.3. Results 

We can now state the result of the computation: 

Proposition. There are 4I smooth lattice polygons with at most 12 lattice points. 



Vertices 


3 


4 


5 


6 7 8 


>9 


Polygons 


3 


30 


3 


401 






The complete list of polytopes can be found in the appendix. 
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3.3. Smooth Polytopes in Dimension Three 
3.3.1. Classifying the Normal Fans 

Minimal Fans 

As in the two-dimensional case, we start with a list of smooth minimal fans from Tadao 
Oda [?]. 

Theorem. (Tadao Oda [7] ) Every smooth fan in dimension three, with at most eight 
rays, can, up to isomorphism, be created by finitely many successive equivariant blowing- 
ups of one of the 19 fans with labels 3^ (S'^A^)', (SM^)", 4^ SM^G^, 3H^5^, 4^5^, 4^62, 
324472^ 33415123^ 32425^62, 3H^5'6\ 3H'5^Q\ {3H'5^'y, {SH^b^^ , (S^S^)', (S^S^)", 
^44g4y ^44g4y/ weights as given in Convex Bodies and Algebraic Geometry. 

Remark. 

• All drawings of such fans in the following chapter are stereographic projections of 
the corresponding triangulation of S"^. 

• The labels give the occurence of vertices with different degrees in the triangulation, 
for example 3^4^6^ means 2 vertices of degree 3, 3 vertices of degree 4 and 2 vertices 
of degree 6. See the pictures later in this chapter. 

Remark. With at most eight rays in the fan, from Eulers formula and the fact the the 
fan is simplicial we know that there are at most twelve full- dimensional cones in the fan 
and thus at most twelve vertices in the polytope. This is the main reason we chose 12 
as upper bound for the number of lattice points in the polytope. 

Even in the two-dimensional case, the blowing-up tree was huge so it comes in handy 
that we can use the two-dimensional classification now to reduce the number of fans 
we have to consider. From the 19 fans of the theorem there will be only five left after 
applying the now following criterion and it will help us a lot later during the computation 
of all blowing-ups. 

Smooth Polygons 

We use the smooth polygons we generated earlier to remove fans whose polytopes would 
have too many lattice points on the boundary. Let V^ik) be the set of all smooth lattice 
fc-gons with at most lattice points. Define the following numbers: 

hik) = mm{\P nZ^\ioT Pe Pjv(A;)} 

iN{k) = min{| int(P) n Z^] for P G VN{k)} 

hN{k) = mm{\P n Z^j - | int(P) n Z^j for P G Pw(fc)} 
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From the results of the previous section we can determine the values shown in the table 
below. 



k 


3 


4 


5 


6 


7 


8 


>9 


li2{k) 


3 


4 


8 


7 


>12 


12 


>12 


ii2{k) 








1 


1 




4 




bi2{k) 


3 


4 


7 


6 




8 





Denote by kr the number of full-dimensional cones in F a ray r is contained in and let 
kp be the maximal kr over all rays r in F. 

All vertices of the polytope defined by the full-dimensional cones around a ray r lie on 
the boundary of the facet defined by the ray r, so for every ray r there is a fc^-gon on 
the boundary of the polytope. Since the polytope is smooth, all these polygons have to 
be smooth, too. This yields the following lower bound for the number of lattice points, 
which depends solely on the combinatorial type of the fan. 

Lemma, (smooth lattice polygon criterion) For a smooth lattice polytope P with 
the normal fan F. 

\PnZ^\>\F^''^\ + {bN{kF)-kF)+ J2 ^^(^-) 

Remark. This bound can be improved by considering the combinatorial type of the fan 
more precisely. However, for our computation it suffices. 

Remaining Minimal Fans 

Now we can check the fans from the theorem for the above criterion, the labels of the 
fans are exactly the numbers k* we are looking for. 

• First, consider all fans with eight rays and twelve full-dimensional cones. That are the 
fans labeled A'6\ 3H^7\ 3H'5'6\ 3H^5^6\ 3H''5'6\ 3H'5^6\ {3H'5'6'y, {3H'5'6'y', 
(3256)', (3256)", (4^5^)' and (4^5^)", which cannot be blown up as this would add two 
additional full-dimensional cones. From the labels one can see that they all have at least 
one ray contained in five, six or seven full- dimensional cones (5*, 6* or 7*), this means 
the corresponding polytope would have a five-, six- or scvcn-gon on the boundary. But 
since any smooth five- or six-gon has at least one interior lattice point and there is no 
smooth seven-gon with less than 13 lattice points at all, the polytope would have at 
least 13 lattice points on the boundary. So none of these fans admits a smooth lattice 
polytope with at most 12 lattice points. 

• Now consider the fans labeled 4^5^ and 3^4^5^ with ten full-dimensional cones. 
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(4^5^) Here we have two non-incident rays each contained in five full- 
dimensional cones, which means we have two edge-disjoint five- 
gons in the polytope. But we know that there is no smooth 
five-gon with less than eight lattice points and thus any smooth 
polytope with this fan as normal fan would have at least 16 ver- 
tices. 

A smooth polytope corresponding to any blowing-up of this fan would have 12 
vertices and either two six-gons or a five-gon and a six-gon and thus have at least 
14 lattice points. 

(3^4^5^) This fan contains three rays with five neighbors, so any polytope 

with this normal fan would have at least three lattice points in au^s^ 
the relative interior of the facets and thus at least 13 lattice 
points in total. Any blowing-up would leave at least one of the 
five-gons unchanged and thus there would still be at least 13 
lattice points, since we have now two additional vertices from 
the blowing-up and one in the interior lattice points in the five- 
gon. 

• By the above criterion, the fans numbered 3"^, (3^4^)', (3^4^)", 4^ and 3^4^6^ could 
appear as normal fans for smooth lattice polytopes with at most 12 lattice points and 
will be the input to the algorithm. 

This leaves us with the following five fans or classes thereof with rays and edge-para- 
meters as drawn [7]. 




oo = (0,0,1) 



-1 

"-^(0,1,0) 

A 

-1 -1 



-1 



-1 



(-1,-1,-1) 
-1^ 



(1,0,0) 

^-1 



(3^43)' oo = (0,1.0) 
1 

\ / 



(0,0,-1)1 " (1,0,0) 
-1 



-1 



(-1,-1, -a) 


/ 







-1 — -a- 



(0,0, 1) 



-1 
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(3^4^)" cx) = (0,0,1) 
-1 c 

\ / 



1- 

{-l,b,c)\ /I (0,1,0) 

b-c -1 



-1 




(0,-1,-1) 
-b 

/ 
c-b 




-1—0 



(1,0,0) 



46 00 = (0,0,1) 

b c 



(-l,-a,c) 



— 



(0,1,0) 



-6 -c 

I 

ab- c -b 
(0,-1,6)1/" ^ _ ^ \| (1,0,0) 



c — ab 







324852 



a + 1 

(0,0,-1) N^* 



-1 a+l 



00 = (1,0,0) 








(0,0,1) 



-a -1 







(0,1,-1) 
2 ' 

/ 

-1 



2 1 



-1—2 



(0,1,0) 
2 



-1 



Tree of Blowing-ups 



Again, the first step of the algorithm is to generate all possible normal fans by successive 
blowing-ups. There are now two possible cases, blowing up in a 3-dimensional cone or 
a 2-dimensional cone, which corresponds in the polytope to cutting off a vertex or an 
edge respectively, as shown in figure El 




Figure 2: Blowing-up of a three-dimensional fan and the cutting-off in the corresponding 
polytope 



Generating all blowing-ups successively gives again a huge rooted tree, but we now have 
more ways to reduce the number of nodes we have to check. 
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(1) We can again limit the depth of the tree since we want at most 12 full-dimensional 
cones and each blowing-up increases the number of such cones by two. 

(2) Blowing-ups in full-dimensional cones can again be done in any order. Let F 
contain k full-dimensional cones. When blowing up in the cone Cj the three new 

cones get the numbers i, A; + 1, A; + 2, so we can ignore all full-dimensional cones 
with index strictly less than i at all nodes in the subtree of this blowing-up. 

(3) For codimcnsion 1 cones there are two cases. Let C be a codimension 1 cone 
contained in two full-dimensional cones Ci and €2- 

(3a) Blowing up 6(F, C) in another codimension 1 cone C", which is contained in 
Ci or C2, will yield a different result than blowing up h{F, C) in the cone C. 
All cones have to be considered in this case. 

(3b) Blowing up h{F, C) in a codimension 1 cone C", which is not contained in Ci 
or C2, will yield the same result as blowing up h{F, C) in C . So we can use 
an ordering on the codimension 1 cones to skip cones with lower index if they 
do not fall into case (3a). 

(4) Blowing up in a full-dimensional cone and later in one of its codimension 1 faces 
yields the same rays and the same combinatorial type as blowing up in the codi- 
mension 1 cone first and then in on new codimension 1 cone as shown below. 

A 7 ''2 f. 7 ''2 




And we have the following new rays: 

d-l d-1 

d-1 d-l 

Xi — ri-\- rii x'2 = x'l + ri — ri + rii — Xi 

i=l 1=1 

So only one of these cases has to be considered. 

(5) As already mentioned, we can use the smooth lattice polygon criterion described 
earlier. Fans that violate this criterion will be considered for further blowing-ups, 
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but will not be kept for the later steps of the algorithm. We cannot cut off the 

tree here because the criterion is not monotone: 

Assume we have a polytope with one vertex that is contained in three five-gons, 
then cutting-off this vertex will replace all with six-gons. So in total this new 
polytope will have two more vertices but may have less lattice points than the 
original polytope, since the smallest smooth five-gon has 8 lattice points in total 
and the smallest smooth six-gon is the hexagon with only 7 total lattice points. 

This has been implemented with one variable for each full-dimensional and each codi- 
mension 1 cone that can have the values 'always ignore', 'ignore', 'consider', 'always 
consider'. We will create the blowing-up and continue the recursion only for cones with 
value 'consider' or 'always consider'. The specific cases are: 

(2) Visited full-dimensional cones are marked 'ignore' for the subtree. 

(3a) Visited codimension 1 cones are marked 'ignore' for the subtree if they haven't 
been marked 'always consider' or 'always ignore'. 

(3b) Codimension 1 cones are marked 'always consider' if they were on 'ignore' before 
and are contained in one of the full-dimensional cones containing the codimension 
1 cone of the current blowing-up. 

(4) Codimension 1 cones are marked 'always ignore' if they were contained in the 
full-dimensional cone of the current blowing-up. 

(x) Newly created cones are always set to 'consider'. 
Pcirameters for the Minimal Fans 

Most of the remaining minimal fans are actually classes of fans with parameters a, 6, c G 
Z. Hence we have to find bounds for their parameters and find ways to use those bounds 
for their blowing-ups. 

We can find bounds for their parameters analogously to the two-dimensional case: When- 
ever a parameter x appears as ac,i — x and acj = —x we know —N < x < N and 
get: 

[3^] has no parameters, 
[(3^4^)'] < a < using the symmetry of the fan. 



[(3^4=^)"] 
[4^] 
[3^4^6^] 



-N <h,c<N, 
- N <a,b,c<N, 
here we use ac,i — ' 



2a + 1 and acj — —2a — 1, so 



-N-l 
2 



<a< 



N-l 
2 
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Pcirameters in the Blowing-up 

In three dimensions we have to consider the following two different cases of blowing-ups: 

• When blowing-up in a full-dimensional cone we can do the same as in the two- 
dimensional case. 

Denote by Ci = {^i, ni, 722}, C2 = {f'2,ni,n2} two full-dimensional cones contain- 
ing a codimension 1 cone C = {ni,n2} with ri + r2 = ac,i^i + 00,2^2- We blow 
up in the cone C2, so the new full-dimensional cones containing C are Ci and 
^2 — {'^2' ^2} with r2 — r2 + rii + n2- We obtain 

ri + r'2 = ac,ini + ac,2?^2 + Ui + n2 = {ac,i + l)ni + {ac,2 + l)ri2 

and this gives ac,i+l > —N and thus ac,i > —N—1. That means we have to loosen 
the bounds on a parameter by one if it appears at one of the codimension 1 cones 
contained in the cone of the current blowing-up, just as in the two-dimensional 
case. 

• For the blowing-up in a codimension 1 cone choose Ci, C2 and C as in the full- 
dimensional case. When blowing up in the cone C at n = ni + n2 we get six 
new cones C[ = {ri,ni,n}, C'2 = {r2,ni,n} containing C = {ni,n} and C" = 
{ri,n2,n}, = {r2,n2,n} containing C" = {n, 112}. Now check the parameters 
at the new cone C: 

ri + r2 = ac,ini + 00,2^2 = ac,ini + ac,2{n - rii) = {ac,i - ac,2)ni + ac,2n 

and at C" we get equivalently: 

ri + r2 = ac,ini + 00,2^2 = ac,i{n - 712) + 00,2^2 = ac,in + {ac,2 - ac,i)^2 

So, in the blowing-up, we have found two codimension 1 cones each having one of 
the edge-parameters of the original fan, ac',2 = 0-0,2 and ac",i = 0,0,1- These edge- 
parameters of the two new cones C", C" yield the same bounds for the parameters 
of the fan as the codimension 1 cone C in the original fan. Which means we can 
just pass on the bounds to the new blowing-up. 

So we have again all neccessary bounds for the parameters of all fans and can create a 
hst of completely defined smooth fans that serve as normal fans for our polytopes. 

3.3.2. Creating the Polytopes 

We can now continue as in the general and the two-dimensional case using polymake 
in combination with 4ti2, normaliz2 and Latte macchiato to compute all smooth 
polytopes with at most 12 lattice points. 
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3.3.3. Results 



Running this algorithm for the remaining five minimal smooth fans yields the following 
result. 

Proposition. There are 33 smooth three-dimensional lattice polytopes with at most 12 
lattice points. 



Vertices 


4 6 8 >10 


Polytopes 


2 25 6 



Again the list of polytopes can he found in the appendix. 



4. Final Comments 

We now have a list of smooth lattice polytopes in dimension two and three with at most 
12 lattice points. For dimension three the bound 12 may seem rather low and in fact 
none of the 33 lattice polytopes from this classification has an interior lattice point. So 
these polytopes may not bring great insight, but this is a first step and there are several 
points in the algorithm where it could be improved. 

By implementing a different way to directly generate all smooth normal fans one could 
skip the big recursion calculating all blowing-ups, which is one of the most time-consu- 
ming steps in the classification algorithm. And one could then overcome the limits of at 
most 12 vertices and thus lattice points and dimension at most three. 

The second point where one could work on is the calculation of lattice points of the 
polytope containing all right-hand sides. The dimension of this polytope is equal to 
the number of rays of the fan minus the dimension. So when we want to reach more 
complex smooth polytopes we will have rather high-dimensional polytopes containing 
the right-hand sides, so we might need better bounds for the right-hand sides and/or 
faster algorithms for computing the lattice points. 
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A. List of Smooth Polygons 




B. List of Smooth Polytopes 




(0,0,0) 
(0,0,1) 
(1,1,0) 



(0,0,0) 
(0,0,1) 
(1,1,0) 



(0,0,0) 
(0,0,1) 
(2,1,0) 



(0,1,0) 
(5,0,0) 
(1,0,1) 




(0,1,0) 
(4,0,0) 
(1,0,1) 



(0,0,0) (0,1,0) 
(0,0,1) (4,0,0) 
(3,1,0) (1,0,1) 



(0,1,0) 
(4,0,0) 
(2,0,1) 



(0,0,0) (0,1,0) 
(0,0,1) (2,0,0) 
(2,1,0) (1,0,1) 



(0,1,0) 
(0,0,0) 
(1,0,1) 



(1,0,0) 
(0,0,1) 
(0,1,1) 



(0,1,0) 
(0,0,1) 
(2,0,0) 
(1,0,1) 



(0,0,0) 
(0,1,1) 
(3,1,0) 
(2,1,1) 



(0,1,0) 
(0,0,1) 
(2,0,0) 
(1,0,1) 




(0,0,0) (0,1,0) 
(0,0,1) (6,0,0) 
(1,1,0) (1,0,1) 




(0,0,0) 
(0,0,1) 
(2,1,0) 



(0,1,0) 
(5,0,0) 
(1,0,1) 




(0,0,0) 
(0,0,1) 
(2,1,0) 



(0,1,0) 
(5,0,0) 
(2,0,1) 



(0,0,0) (0,1,0) 
(0,0,1) (4,0,0) 
(4,1,0) (1,0,1) 




(0,0,0) 
(0,0,1) 
(2,1,0) 



(0,1,0) 
(3,0,0) 
(1,0,1) 



(0,0,0) 
(0,0,1) 
(3,1,0) 



(0,1,0) 
(4,0,0) 
(2,0,1) 




(0,0,0) 
(0,0,1) 
(3,1,0) 



(0,0,0) 
(0,0,1) 
(2,1,0) 



(0,0,0) 
(0,0,1) 
(1,1,0) 



(0,1,0) 
(1,0,0) 
(1,0,1) 

(1,1,1) 



(0,1,0) 
(0,0,1) 
(2,0,0) 
(2,0,1) 



(0,1,0) 
(5,0,0) 
(1,0,1) 




(0,1,0) 
(4,0,0) 
(1,0,1) 




(0,1,0) 
(3,0,0) 
(1,0,1) 




(0,1,0) 
(0,1,1) 
(2,0,0) 



(0,0,0) (0,0,1) 
(2,0,0) (0,2,0) 
(1,0,1) (0,1,1) 



(0,0,0) 
(0,1,1) 
(2,1,0) 
(2,1,1) 



(0,0,0) (0,1,0) 
(0,0,1) (3,0,0) 
(3,1,0) (1,0,1) 




(2,0,0) (0,2,0) 
(2,0,1) (0,2,1) 
(0,3,0) (3,0,0) 



(0,1,0) 
(0,0,1) 
(3,0,0) 
(1,0,1) 



(0.0.0) 
(0.1.1) 
(3,1,0) 
(1,1,1) 
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C. polymake extension 

C.l. Parameterized Fans (pf an. rules) 



# a parameterized class of fans 
declare object ParameterFan <Ratlonal > { 

# maximal number of lattice points 
property POINT.LIMIT : Int ; 

# dimension of the fan 
property DIM : Int ; 

# incidence sets (of the rays) defining the full -dimensional cones or 

# vertices of the polytope 
property VERTICES : Array<Set>; 

# number of vertices 
property N.VERTICES : Int; 

# incidence sets (of the vertices) defining the codimension 1 cones or 

# edges between vertices of the poly topes 
property EDGES : Array<Set>; 

# number of edges 
property H.EDGES : Int; 

# rays of the fan 

property RAYS : Matrix <Rational > ; 

# number of rays 
property K.RAYS : Int; 

# number of different parameters 
property N.PARAMETER : Int ; 

# rows with positions in RAYS matrix, parameter number and coefficient 
property PARAMETER.POS : Matrix<Int>; 

# lower and upper bound for each parameter (rowwise) 
property PARAMETER.BOUNDS : Matrix<Int>; 

# searches full - dimensional cones sharing a common codimension 1 cone 
rule EDGES : VERTICES , N.VERTICES , DIM { 

my $vert = $this ->VERTICES ; 
my Sedges = () ; 
local $_ ; 

# check all possible pairs 

foreach ( all.subset s_of _k (2 , . . ( $this ->N_VERTICES -1) ) ) { 
my $r = $vert ->[$_-> [0] ] * $vert ->[$_-> [1] ] ; 
if ($#$r == ($this->DIM - 2)) { 
push Sedges, new Set($_); 

} 

} 

$this->EDGES = new Array <Set >( Sedges ) ; 

} 

rule DIM : RAYS { 

$this->DIM = $this ->RAYS ->cols ; 

} 

rule N_EDGES : EDGES { 

$this ->N_EDGES = $#{ $this ->EDGES } + 1; 
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} 



rule N_VERTICES : VERTICES { 

$this ->N_VERTICES = $#{ $this -> VERTICES } + 1; 

} 

rule N_RAYS : RAYS { 

$this ->N_RAYS = $this ->RAYS ->rows ; 

} 

# main method creating the list of polytopes for this fan 

# calls all necessary other functions 

# returns: the polytopes as perl array of Polytope objects 

# input: integer value determining the print- or debug-level 
user.method POLYTOPES { 

my ($thls , $print) = S. ; 
my Opoly = () ; 

print "Creating blowing-ups ...\n" if($print); 
my Sblowups = $this -> create _blowups ( Sprint ) ; 
print "Done creating blowing -ups . \n " if($print); 
my $i=l; 

print "Generating polytopes ...\n" if (Sprint); 
foreach my $pfan (Oblowups) { 

print "blowup " , $i , " of ", s calar ( Sblowups )," \n" if (Sprint); 

my afans = $pf an -> get f ans ( ) ; 
my $j = 1; 

foreach my $fan (Sfans) { 

print "\tfan " , $i , " : " , $ j , " of ", scalar ( Oblowups ),":" , 

scalar Of ans ), "\n" if (Sprint); 
my $num = $f an->N_POLYTGPES ; 
if ($num > 0) { 

push apoly, $fan->POLYTOPE_LIST ; 

} 

$j++; 

} 

$i++; 

} 

print "Done generating polytopes . \n" if($print); 

print "Found " , s calar ( Spoly ) , " polytopes including duplicates . \n" if($print); 
print "Removing duplicates ...\n" if($print); 
my @final_poly = (); 

outer: for (my $i = s calar ( Opoly ) - 1 ; $i >= 0; $i--) { 
my $p = Spoly [$i] ; 

for (my $j = $i-l; $j >= 0; { 
my $q = Spoly [$j] ; 

# reduce neccessary isomorphy checks by checking number of vertices first 
if ($p->N_VERTICES == $q->M_VERTICES ) { 

if ( lattice_isomorphic_smooth_polytopes ($p , $q) ) { 
next outer ; 

} 

} 

} 

push 9f inal_poly , $p ; 

} 

print "Found ", scalar (Sfinal_poly) , " non - i somorphic polytopes with this normal fan.\n"; 
return Sfinal_poly; 

} 

# starts the recursion generating all blowing-ups returning them in an perl array 

# printlimit specifies up to which depth to print messages 
user_method create_blowups { 

my ( $f an , Sprint 1 imlt ) = S_; 

# incidence vectors storing which vertices and edges should be blown-up in the tree, 

# zero vectors at the beginning 

my $vertexlist = new Vector < Int >( $f an ->POINT_LIMIT ) ; 
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my $edgelist = undef ; 

$edgelist = new Vect or < Int > ( ( $f an - > DIM * $f an - > PO INT.LIMIT ) /2) if ($fan->DIM > 2); 



my (Sblowups = r e cur s ion ( $f an , $vert exl ist , $edgelist , , $pr int limit ) ; 
return Sblowups ; 

} 

# creates a new fan for every combination of parameter values and 

# returns the list of fans as perl array 
user.method getfans { 

my $bf = shift; 

my $rays = $bf->RAYS; 

my af ans = () ; 

if ($bf ->N_PARAMETER > 0) { 

my $npar = $bf ->N_PARAMETER ; 

my $pb = $bf ->PARAMETER_BOUNDS ; 

# vector storing the current parameter values 

my $par = new Vector <Int >( $pb -> col (0) ) ; 

my $cont = 1 ; 

do { 

my $r = new Matrix <Rational >( $rays ) ; 

# add the parameter values to every position it occurs 
foreach my $pos ( (a{$bf - > P ARAMETER.POS }) { 

$r->($pos -> [0] , $pos -> [1] ) += $par -> [$pos -> [2] ] * $pos->[3]; 

} 

push ©fans ,new Fan < Rat i onal > ( POINT_LIMIT = > ( $bf - > PO INT_LIMIT ) , DIM=>( $bf ->DIM) , 
RAYS=>$r , VERTICES =>($bf-> VERTICES) , EDGES =>( $bf ->EDGES ) ) ; 



# generate next vector of parameter values 

# increase first parameter and check bounds , ■ . . 
$par->[0] ++; 

for (my $i =0; $i < $npar ; $i++) { 
if ($par->[$i] > $pb->[$i]->[l]) { 
if($i == $npar-l) { 

$cont = 0; 
} else { 

$par->[$i] = $pb->[$i] ->[0] ; 
$par->[$i+l] ++; 

} 

} 

} 

} while ($cont == 1) ; 
y else -[ 

push (Sfans ,new Fan <Rat ional > ( POINT.LIMIT => ( $bf -> POINT.LIMIT ) , DIM=>( $bf ->DIM) , 
RAYS=>$rays , VERTICES =>( $bf -> VERTICES ) , EDGES =>( $bf ->EDGES ) ) ; 

} 

return Sfans ; 

} 



# main recursion creating the blowing-ups 
sub recursion { 

my ( $f an , $vertexlist , $edgelist , $depth , $pr intlimit ) = a_ ; 

my @b = (); 

if ($fan->DIM > 2) { 

# in dimension > 2 check smooth polygon criterion 
if (checkpolygons ($f an) <= $f an ->P0 INT.LIMIT ) { 
push ab , $fan; 

} 

} else { 

push @b , $f an ; 

} 



# if blowing-up does not exceed maximal number of vertices 
if ($f an->N_VERTICES <= $f an->POINT_LIMIT - ( $f an ->DIM - 1 ) ) { 
# blowing-up of full -dimensional cones 
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for(my $i = 0; $i < $f an -> N_ VERTICES ; { 
if ($vertexlist -> [$i] < 1) { 

print "do ver t exblowup : " , $depth , " - " , $i , " \n " if ($depth < $pr intlimit ) ; 
my ( $newf an , $elist ) = do_ vert ex_blowup ( $f an , $i , $edgel ist ) ; 
push Sb , recur sion ( $newf an ,( new Vect or < Int > ( $ vert exli st ) ) , 

$elist , $depth+l , Sprint limit ) ; 

$ vert exli St -> [$i] =1 ; 

} 

} 

# if dim > 2 blowing-up of codim 1 cones 
if($fan->DIM > 2) { 

for (my $i =0; $i < $f an ->N_EDGES ; $i++) { 
if (Sedgelist ->[$i] < 1) { 

print "do edgeblowup : " , $depth , " -" , $i , " \n" if ($depth < $printlimit ) ; 
my ( $newf an , $vlist , $eli St ) = do_edge_blowup ($f an , $i , $vertexlist , Sedgelist ) ; 
push 9b, recur s i on ( $newf an , $vlist , $elist , $depth+l , Sprint 1 imit ) ; 
if ($edgelist -> [$i] == 0) { 
Sedgelist -> [$i] = 1; 

} 

} 

} 

} 

} 

return 9b; 

} 

# create the new fan after the blowing -up of full-dim cone $i 
sub do_vertex_blowup { 

my ($f an , $i , Sedgelist) = «_ ; 

# set of rays defining the vertex 
my $list = $f an->VERTICES -> [$i] ; 
my $nrays = $f an->N_RAYS ; 

my $elist = undef ; 

$elist = new Vector <Int >( $edgelist ) if ( def ined ( $edgelist ) ) ; 

# create and add the new ray 

my $ray = new Vector <Rational >( $fan->DIM ) ; 

local $_; 

foreach (OSlist) { 

$ray += $f an ->RAYS ->[$_] ; 

} 

my $newrays = $fan->RAYS / $ray ; 

# prepare edges array, vertices array and array with positions for new vertices 
my Snewedges = a{$f an->EDGES}; 

my @pos = ( $i ) ; 

my Svertices = ® { $f an -> VERTICES } ; 

for (my $k =0; $k < $fan->DIM - 1; $k++) { 

push Spos , ($f an->N_VERTICES+$k) ; 

push Svertices, new SetO; 

} 

local $_; 

foreach ( all.subset s_of _k (2 , 9pos ) ) { 
push Snewedges , (new Set($_)); 

} 

my $k = 0; 

my $newvertices = new Arr ay < Set > ( @ vert i ces ) ; 
my $newedges = new Array <Set >( Snewedges ) ; 

# find neighbors of vertex we are cutting off, correct the edges 

# and create new vertices 

for (my $j = 0; $j < $f an - > N_EDGES ; $j+ + ) { 
my $edge = $f an ->EDGES - > [ $ j ] ; 
if ($edge->contains ($i)) { 
my $v = ($edge - $i)->[0]; 
# new vertex 
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$newvertices -> [Opos [$k] ] = ( $f an -> VERTICES ->[ $v] * $list) + $nrays ; 

# set new edge endpoints 

$newedges -> [$j ] = new Set ($v , Spos [$k] ) ; 

# set neighboring edges to always ignore, 

# we do vertex -> neighboring edge via edge -> edge 
$elist->[$j] = 2 if (defined ( $edgelist )) ; 

$k++; 

} 

} 

my $newfan; 

# odd parameters for the new ray and loosen bounds if neccessary 
my $par = new Matrix <Int >() ; 
if ($f an->N_PARAMETER > 0) { 

my Sparam = 8{$f an ->PARAMETER_POS } ; 

local $_; 

my $p_indices = new Vector < Int >( $f an ->N_PARAMETER ) ; 

# find parameters appearing in the rays defining the currently blown -up cone 
foreach ( a{ $f an -> P ARAMETER_POS }) { 

if ($list ->contains ($_ - > [0] ) ) { 

push Sparam , (new Vect or < Int > ( $nr ay s , $_ - > [1] , $_ - > [2] , $_ - > [3] ) ) ; 
$p_ indices ->[$_->[2]]=l; 

} 

} 

$par = new Matr ix < Int > ( Opar am ) ; 

# loosen hounds on all parameters found earlier 

my $pbounds = new Matr ix < Int > ( $f an -> P ARAMETER.BOUMDS ) ; 
for (my $j = 0; $j < $f an - > N_P ARAMETER ; $j++) { 
if ($p_indices -> [$j ] > 0) { 

$pbounds -> [$j ] -> [0] --; 

$pbounds -> [$j ] -> [1] ++; 

} 

} 

# create new fan with all modified data 

$newfan = new ParameterFan<Rational >(POINT_LIMIT=>$f an->PQINT_LIMIT , 

DIM=>$f an->DIM , VERTICES=>$newvertices , EDGES=>$newedges , RAYS=>$newrays , 
PARAMETER_POS=>$par , PARAMETER_BOUHDS=>$pbounds , N_PARAMETER=> $f an ->N_PARAMETER ) ; 

} else { 

$newfan = new Par amet erFan < Rat i onal > ( POINT_LIMIT = > $f an - > PO INT_LIMIT , 

DIM=>$f an->DIM , VERTICES=>$newvertices , EDGES => $newedges , RAYS=>$newrays , 
PARAMETER_POS=>$par , N_PARAMETER=>$f an->N_PARAMETER) ; 

} 

return ( $newf an , $elist ) ; 



# create the new fan after the blowing-up of full-dim cone $i 
sub do_edge_blowup { 

my ( $f an , $i , $vertexlist , $edgelist ) = ®_ ; 

my $vlist = new Vector <Int >( $vertexlist ) ; 

my $elist = new Vector < Int >( $edgelist ) ; 

# indices of the two vertices of the edge 
my $edge = $f an->EDGES -> [$i] ; 

my $el = $edge->[0]; 
my $e2 = $edge->[l]; 

# sets of the rays defining these vertices 
my $vl = $f an->VERTlCES -> [$el] ; 

my $v2 = $f an->VERTlCES -> [$e2] ; 

# two rays defining this edge (as set and separate) 
my $er = $vl*$v2; 

my $erl = $er-> [0] ; 
my $er2 = $er->[l] ; 

# total number of vertices and rays 
my $nv = $f an->N_VERTICES ; 

my $nrays = $f an->H_RAYS ; 

# add new ray 
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my $newrays = $fan->RAYS / ( $f an ->RAYS -> [$er - > [0] ] + $f an ->RAYS - > [$er - > [1] ] ) ; 

# new edg earray 

my Snewedges = (a{$f an -> EDGES } ; 

# create edges to and between the two vertices that will be put at the end 
push Snewedges, new Set ($nv , $nv+l) ; 

push Snewedges, new Set ($el , $nv) ; 
push Snewedges, new Set ($e2 , $nv+l) ; 
my $k = 0; 

# create new vertices 

my $newvertices = new Array <Set >( S{$f an -> VERTICES }, new Set () , new SetO); 
$newvertices -> [$el] = $vl - $er2 + $nrays ; 
$newvertices -> [$e2] = $v2 - $er2 + $nrays ; 
$newvertices -> [$nv] = $vl - $erl + $nrays ; 
$newvertices -> [$nv+l] = $v2 - $erl + $nrays ; 

# set the new vertices to 'consider' for later blowing-ups 
$vlist -> [$el] = 0; 

$vlist -> [$e2] = 0; 

# create the new edges and update all old edges 

# we skip all just added edges and the cut-off edge 
my $newedges = new Array <Set >( Snewedges ) ; 

for (my $j = 0; $j < $f an - > N_EDGES ; $j++) { 
next if ($i==$j ) ; 
my $e = $f an->EDGES -> [$ j ] ; 

if ($e->contains ($el) I I $e -> contains ( $e2 ) ) { 
my $v = ($e - $edge)->[0]; 

if (! ($fan->VERTICES->[$v]->contains($erl))) { 

Snewedges ->[$j] = new Set($v,$nv) if ($e -> contains ( $el )) ; 
Snewedges ->[$j] = new Set($v,$nv+1) if ($e ->contains ( $e2 ) ) ; 

} 

$elist->[$j] = -1 if ($elist -> [$j ] ==1) ; 

} 

} 

my $newf an ; 

# check if the new ray has parameters 
my $par = new Matr ix < Int > () ; 

if ($f an->N_PARAMETER > 0) { 

my Sparam = ®{$f an ->PARAMETER_PDS } ; 
local $_; 

# search for parameters for the rays defining the new one 
foreach ( a{ $f an -> P ARAMETER_POS }) { 
if ($er->contains ($_ -> [0] ) ) { 

push Sparam , (new Vect or < Int > ( $nrays , $_ - > [1] , $_ - > [2] , $_ - > [3] ) ) ; 

} 

} 

$par = new Matrix <Int >( Sparam) ; 

$newfan = new ParameterFan<Rational >(POINT_LIMIT=>$f an->POINT_LIMIT , 

DIM=>$f an->DIM , VERTICES=>$newvertices , EDGES=>$newedges , RAYS=>$newrays , 
PARAMETER_POS=>$par , PARAMETER.BOUNDS =>$f an ->PARAMETER_BOUNDS , 
N_PARAMETER=>$f an->N_PARAMETER) ; 

} else { 

$newfan = new Par amet erFan < Rat i onal > ( POINT_LIMIT = > $f an - > PD INT.LIMIT , 

DIM=>$f an->DIM , VERTICES=>$newvertices , EDGES=>$newedges , RAYS=>$newrays , 
N_PARAMETER=>$f an->N_PARAMETER) ; 

} 

return ( $newf an , $ vlist , $elist ) ; 



# check the smooth polygon lower bound for the number of lattice points 
sub Che ckpolygons { 

my ($fan) = S_ ; 

my $rays = $f an->N_RAYS ; 
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# results from the 2- dimensional classification ^ 

# "13" means not possible for <= 12 lattice points 

# FIXME for more than POINT_LIMIT > 12 

my $innercounts = new Vector <Rational >(' 1 1 13 4 13 13'); 
my $bordercounts = new Vector <Rat ional >(' 3 4 7 6 13 8 13 13'); 

# collect degrees of all rays (number of vertices they appear in) 
my $sizes = new Vector<Int>( '0 ' x $rays ); 

local $_; 

foreach (a{$f an->VERTICES}) { 
foreach my $x (a$_) { 
$sizes->[$x] ++; 

} 

} 

# find maximal degree, and sum up all interior lattice points 
my $count = 0; 

my $inax = 0; 
local $_; 

foreach ((S$sizes) { 
if ($_ > $max) { 
$max = $_ ; 

} 

next if ($_ > 8) ; 

$count += $innercounts -> [$_ -3] ; 

} 

if ($max > 8) { 
return 13; 

} 

# odd boundary lattice points for the largest polygon 
$count += $bordercounts -> [$max -3] ; 

$count += $f an->N_VERTICES - $max ; 
return $count ; 



C.2. Fans (fan. rules) 



# a completely defined fan 
declare object Fan<Rational > { 

# maximal number of lattice points 
property POINT.LIMIT : Rational; 

# dimension of the fan 
property DIM : Int ; 

# incidence sets (of the rays) defining the full - dimensi onal cones or 

# vertices of the poly tope 
property VERTICES : Array<Set>; 

# number of vertices 
property N_VERTICES : Int; 

# incidence sets (of the vertices ) defining the codimension 1 cones or 

# edges between vertices of the poly topes 
property EDGES : Array<Set>; 

# number of edges 
property K.EDGES : Int; 

# rays of the fan 
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property RAYS : Matrix <Rational > ; 

# number of rays 
property N_RAYS : Int ; 

# polytope approximating the set of valid right-hand sides 
property RHS.POLYTOPE : Polytope <Rational > ; 

# set of right-hand side vectors defining polytopes 

# with at most POINT_LIMIT lattice points 
property RHS.VECTORS : Set < Vector < Integer >> ; 

# set of right-hand side vectors defining polytopes 

# with at most POINT_LIMIT lattice points 
property N.POLYTOPES : Integer; 

# returns a perl array of all polytopes with this normal fan and 

# at most POINT_LIMIT lattice points 
user_method POLYTOPE_LIST { 

my $fan = shift; 
my (Slist = (); 

# create polytope for every rhs vector 
foreach my $rhs ( a{ $f an -> RHS.VECTORS } ) { 

my $ineq = (new Vector <Rat ional >( $rhs ) ) I ( $f an ->RAYS ) ; 

my $p = new Polytope <Rational >( INEQUALITIES=>$ineq) ; 

#my $v = $p->N_VERTICES ; 

#my $lp = $p->N_LATTICE_POINTS; 

push Slist , $p ; 

} 

return Slist ; 

} 

# searches full -dimensional cones sharing a common codimension 1 cone 
rule EDGES : VERTICES , N.VERTICES , DIM { 

my $vert = $this ->VERTICES ; 
my Sedges = () ; 
local $_; 

# check all possible pairs 

foreach ( all _subs et s _ of _k (2 , ..( $this -> N.VERTICES - 1) ) ) { 
my $r = $ vert ->[$_-> [0] ] * $vert ->[$_-> [1] ] ; 
if ($#$r == ($thls->DIM - 2)) { 
push Sedges, new Set($_); 

} 

} 

$this->EDGES = new Array <Set >( Sedges ) ; 

} 

rule DIM : RAYS { 

$this->DIM = $this->RAYS->cols ; 

} 

rule N.EDGES : EDGES { 

$this ->N_EDGES = $#{$this ->EDGES} + 1; 

} 

rule N.VERTICES : VERTICES { 

$this ->N_VERTICES = $#{ $thi s -> VERTICES } + 1; 

} 

rule N_RAYS : RAYS { 

$this ->N_RAYS = $this ->RAYS ->rows ; 

} 

rule N.POLYTOPES : RHS.VECTORS { 

$this ->N_POLYTOPES = $#{ $this ->RHS_VECTORS }+l ; 

} 
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# searches in RHS_POLYTOPE for RHS-vectors defining poly topes 

# with at most POINT_LIMIT lattice points 

rule RHS_VECTQRS : RHS_POLYTOPE , RAYS , POINT.LIMIT , DIM { 
my $valid = $this ->RHS_POLYTDPE ; 
my $limit = $this ->POINT_LIMIT ; 
my $dim = $this->DIM; 
my $rays = $this->RAYS; 

# check the polytope for feasibility and number of lattice points first 
if ($valid->FEASIBLE == 0) { 

$this ->RHS_VECTORS=new Set < Vector < Integer >>( ) ; 
return ; 

} 

if ($valid->N_LATTICE_POINTS == 0) { 

$this ->RHS_VECTORS=new Set < Vector < Integer >> () ; 
return ; 

} 

# remove first coordinate (homogenized) 

my $rhs = $ valid - >LATTICE_POINTS - >minor (All , range ( 1 , $r ays ->rows ()- $dim )) ; 
my arhs = a$rhs ; 
Srhs = sort (Srhs ; 

my @good_rhs = () ; 
my $i=0; 
my $sched2 ; 

while (Srhs > 0) { 

# create new RHS vector and inequality matrix 

my $x = ((zero.vector <Int>($dim)) I shif t ( 9rhs ) ) ; 
my $ineq = (new Vector <Rational >( $x ) ) I $rays ; 

# define polytope and check number of lattice points 
my $p = new Polytope <Rat ional >( INEQUALITIES => $ineq) ; 

if ($p->N_LATTICE_PDINTS <= $limit) { 

# found a new polytope 
push Sgood.rhs , $x; 
$i++; 

} else { 

# we only keep RHS-vectors with at least one entry with smaller 

# absolute value than the current vector 
my anewrhs = () ; 

point: for my $v (Srhs) { 

for(my $j = 0; $j < $rays -> r ows () - $dim ; $j+ + ) { 
if (abs ($v-> [$j ] ) < abs($x->[$j+$dim])) { 
push 8newrhs,$v; 
next point ; 

} 

} 

} 

arhs = Qnewrhs ; 

} 

} 

$this ->RHS_ VECTORS =new Set < Vector < Integer >>( Sgood.rhs ) ; 

} 

# creates the polytope approximating the set of right-hand sides 

rule RHS_POLYTOPE : RAYS , VERTICES , EDGES , POINT.LIMIT , N.EDGES , N.VERTICES { 
my $vert = $this ->VERTICES ; 
my Sedges = $this ->EDGES ; 
my $rays = $this->RAYS; 
my $dim = $r ays -> cols () ; 
my $r = $rays ->rows () ; 
my $limit = $this ->POINT_LIMIT ; 
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# stores the sum of all edge - lengths 

my $totallength = new Vector <Rat ional >( $r+l) ; 

# first stores 1(e) >= 1 inequalities , second one inequalities for l(e)<=... 
my aineq = () ; 

my ai.ineq = () ; 

# calculate all edge lengths and create inequalities 
foreach (@$edges) { 

# vertices defining the edge 

my $fromvert = $vert ->[$_-> [0] ] ; 
my $tovert = $vert ->[$_-> [1] ] ; 

# rays defining the edge (n_l , . . . , n_{d-l}) 
my $neighbors = $f romvert * $tovert ; 

# rays limiting this edge (r_l , r_2) 
my $from = $f romvert - $neighbors; 
my $to = $tovert - $nelghbors ; 

$f rom = $from-> [0] ; 
$to = $to -> [0] ; 

# calculate inequalities 

my ( $edgelimit , $pointlimit ) = edgel imit s ( $ray s , $f rom , $to , $limit , $ne ighbor s ) ; 

$totallength += $pointlimit ; 
push Qineq , $edgelimit ; 
push ai.ineq, $pointlimit ; 

} 

# right-hand side for total length inequality 

$totallength -> [0] = $limit + $this ->N_EDGES - $this ->N_VERTICES ; 
push ai_ineq , $totallength ; 

my $ineq = new Matrix <Rat ional > (\ fiineq ) ; 

my $l_ineq = new Matrix <Rational > (\91_ineq) ; 

# for all rays with all entries positive we know that the right-hand side must 

# be postive and vice versa (we want the negative RHS -vectors ) 
my $signs = new Matr ix <Rat ional >( $r , $r+l ) ; 

my Osigns = () ; 

row: for (my $i =0; $i < $r ; $i++) { 

my $sign = 0; 

for (my $j = 0; $j < $dim ; $ j ++) { 
if ($rays ->($i , $j ) < 0) { 
next row if ($sign > 0); 
$sign = -1; 
} elsif ($rays ->($i , $j ) > 0) { 
next row if ($sign < 0); 
$sign = 1; 

} 

} 

$signs ->($i , $i+l) = -$sign; 



# put all inequalities together and remove the first d columns after the rhs entry 

# because those variables are always zero (we do not want translations) 
$ineq = $ineq / $signs; 

$ineq = $ineq - >minor ( All , range (0 , )+ range ( $dim+ 1 , $r )) ; 
$l_ineq = $l_ineq ->minor ( All , range (0 , 0) + range ( $dim+l , $r )) ; 

my $limited = (new Matrix <Rational >( $ineq) ) / (new Matr ix <Rat ional >( $l_ineq) ) ; 
$this ->RHS_POLYTOPE = new Polytope <Rat ional > ( INEQUALITIES=>$limited ); 
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# calculates the inequalities in terms of the rhs vector given by l<=l(e)<=... 
sub edgelimlts { 

my ($rays , $f rom , $to , $limit , $nb) = @_ ; 

my $dim = $rays - > cols () ; 

# limiting rays 

my $rl = $rays ->roM ( $to ) ; 
my $r2 = $rays ->row ( $f rom ) ; 

# rays defining the edge 

my $N = $rays ->minor ( $nb , All ) ; 

# determine edge-parameters a_{C,i} 

$H = new Matrix <Rational > (transpose ($N) ) ; 
my $v = new Vector <Rational >( $rl+$r2 ) ; 
my $a = lin.solve ($H , $v) ; 
my $sum = 0; 

# and sum them up for the right-hand side 
$sum += $a->[$_] for (0 .. ($dim-2)); 

# vectors storing the edg el eng th >= 1 and <= . . . inequalities 
my $edgelimit = new Vect or <Rat ional > ( $ray s -> r ows ( ) + 1); 

my $pointlimit = new Vect or <Rat ional >( $rays -> rows () + 1); 

# rhs of the inequalities 
$edgelimit -> [0] = (-1); 

$pointlimit -> [0] = (($limit - $sum ) / $dim - 1 ) ; 

# set coefficients 

for (my $i = 0; $i < $dim-l; $i++) { 
my $k = $nb->[$i] + 1; 
$edgelimit -> [$k] = (-1) * $a->[$i]; 
$pointlimit -> [$k] = $a->[$i]; 

} 

$edgelimit -> [$to+l] = 1; 
$pointlimit -> [$to + l] = -1; 
$edgelimit -> [$f rom+1] = 1; 
$pointlimit -> [$f rom+1] = -1; 

return ( $edgel imit , $po int limit ) ; 

} 



C.3. Additional Script (checkiso) 

# removes all lattice isomorphic poly topes from a list 

# only the first one of an isomorphy class is kept 
sub checkiso { 

my (Olist) = a_ ; 
my @res = () ; 

outer: for (my $i = scalar ( Slist ) -1 ; $i >= 0; $i--) { 

my $p = aiist [$i] ; 

for (my $j = $i-l; $j >= 0; $j — ) { 
my $q = aiist [$j ] ; 

if ($p->N_VERTICES == $q->N_VERTICES ) { 

if ($p->N_LATTlCE_P01NTS == $q->N_LATTlCE_P01NTS ) { 

next outer if ( latt i ce_isomorphic_smooth_polyt opes ( $p , $q) ) ; 

} 

} 

} 

push 9res , $p ; 

} 

return iSres ; 
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C.4. Minimal Fans as polymake files 



C.4.1. Dimension Two 



Projective Space 



<?xml version=" 1 . 0" encoding=" utf -8" ? > 



<object naine="P" type =" po lyt ope :: Par amet erFan&lt ; Rat i onal &gt ; " version= 
tm="4b0d8c7b" xmlns="http:// www .math.tu-berlin.de/polymake/#3"> 
<property name=" POINT_LIMIT " value="12" /> 
<property name="DIM" value="2" /> 
<property name =" VERTICES " > 
<in> 

<v>0 l</v> 
<v>l 2</v> 
<v>0 2</v> 
</m> 
</property > 

<property name="RAYS"> 
<m> 

<v>l 0</v> 
<v>0 l</v> 
<v>-l -l</v> 
</m> 
</property > 

<property name = " N_RA YS " value="3" /> 
<property name=" N_VERTICES " value="3" /> 
<property name = " N_P ARAMETER " value="0" /> 
</ obj ect > 



"2.9.7" 



Hirzebruch surfaces 

<?xml version=" 1 . 0" encoding=" utf -8" ? > 

<object name="H" type=" polytope : : ParameterFan&lt ; Rational&gt ; " ver s ion= " 2 . 9 . 7 " 
tm="4af96f7b" xmlns="http:// www .math.tu-berlin.de/polymake/#3"> 
<property name = " PO INT_LIMIT " value="12" /> 
<property name="DIM" value="2" /> 
<property name=" VERTICES " > 
<m> 

<v>0 l</v> 
<v>0 2</v> 
<v>l 3</v> 
<v>2 3</v> 
</m> 
</property > 

<property name="RAYS"> 
<m> 

<v>l 0</v> 
<v>0 l</v> 
<v>0 -l</v> 
<v>-l 0</v> 
</m> 
</ property > 

<property name = " N_RAYS " value="4" /> 
<property name = " N_VERTICES " value = "4" /> 
<property name=" PARAMETER.POS "> 
<m> 

<v>3 1 l</v> 
</m> 
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</ property > 

<property name = " PARAMETER.BOUNDS " > 
<in> 

<v>0 12</v> 
</m> 
</property > 

<property name =" N.PARAMETER " value="l" /> 
</ obj ect > 



C.4.2. Dimension Three 



Fan 3^ 

<?xml version=" 1 . 0" encoding=" utf -8" ? > 

<object name="f" type=" polytope : : ParameterFan&lt ; Rational&gt ; " ver s ion = " 2 . 9 . 7 " 
tm = "4b0d8cld" xiiilns = "http:/ / www .math.tu-berlin. de /polymake /#3 " > 

<property name=" POINT_LIMIT " value="12" /> 

<property name="DIM" value="3" /> 

<property name=" VERTICES " > 
<m> 

<v>0 1 2</v> 
<v>0 1 3</v> 
<v>0 2 3</v> 
<v>l 2 3</v> 
</m> 
</property > 

<property name="RAYS"> 
<m> 

<v>l 0</v> 
<v>0 1 0</v> 
<v>0 l</v> 
<v>-l -1 -l</v> 

</m> 
</prop6rty > 

<property name = " N_RAYS " value="4" /> 
<property name=" N.VERTICES " value="4" /> 
<property name =" N.PARAMETER " value="0" /> 
</ obj ect > 



Fan (3^42)' 

<?xml version=" 1 . 0" encoding= " utf -8 " ? > 

<object name="f" type=" polytope : : ParameterFan&lt ; Rational&gt ; " ver s ion = " 2 . 9 . 7 " 
tm="4b0d571c " xmlns=" http : //www . math . tu-berlin . de /polymake /#3 " > 

<property name=" POINT_LIMIT " value="12" /> 

<property name="DIM" value="3" /> 

<property name=" VERTICES " > 
<m> 



<v>0 


1 


2</v> 


<v>0 


1 


3</v> 


<v>0 


2 


4</v> 


<v>0 


3 


4</v> 


<V>1 


2 


4</v> 


<V>1 


3 


4</v> 



</m> 



</ property > 

<property name="RAYS"> 
<m> 

<v>l 0</v> 
<v>0 1 0</v> 
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<v>0 l</v> 
<v>0 -l</v> 
<v>-l -1 0</v> 
</m> 
</property > 

<property name = " N_RAYS " value="5" /> 
<property name=" N_VERTICES " value="6" /> 
<property name=" PARAMETER.POS "> 
<m> 

<v>4 2 -l</v> 
</m> 
</property > 

<property name=" PARAMETER.BOUNDS " > 
<m> 

<v>0 12</v> 

</m> 
</ property > 

<property name = " N.PARAMETER " value = "l" /> 

</ obj ect > 

Fan (3242)" 

<?xml version=" 1 . 0" encoding= " utf -8 " ? > 

<object name="f" type=" polytope : : ParameterFan&lt ; Rational&gt ; " ver s ion = " 2 . 9 . 7 " 
tm="4b0d8bd2 " xmlns=" http : //www . math . tu-berlin . de /polymake /#3 " > 

<property name=" POINT_LIMIT " value="12" /> 

<property name="DIM" value="3" /> 

<property name=" VERTICES " > 
<m> 



<v>0 


1 


2</v> 


<v>0 


1 


3</v> 


<v>0 


2 


3</v> 


<V>1 


2 


4</v> 


<V>1 


3 


4</v> 


<v>2 


3 


4</v> 



</m> 



</ property > 

<property name="RAYS"> 
<m> 

<v>l 0</v> 
<v>0 1 0</v> 
<v>0 l</v> 
<v>0 -1 -l</v> 
<v>-l 0</v> 
</m> 
</property > 

<property name=" M.RAYS " value="5" /> 
<property name=" N.VERTICES " value="6" /> 
<property name=" PARAMETER.POS "> 
<m> 

<v>4 1 l</v> 
<v>4 2 1 l</v> 

</m> 
</property > 

<property name=" PARAMETER.BOUNDS " > 
<m> 

<v>-12 12</v> 
<v>-12 12</v> 
</m> 
</property > 

<property name =" N.PARAMETER " value="2" /> 
</ obj ect > 
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Fan 4^ 



<?xml version=" 1 . 0" encoding=" utf -8" ? > 

<object name="f" type=" polytope : : ParameterFan&lt ; Rational&gt ; " version=" 2 . 9 . 7" 
tm="4b0d8a3c " xmlns="http : //www . math . tu-berlin . de /polymake /#3 " > 
<property name=" PDINT.LIMIT " value="12" /> 

<property nanie = "DIM" value = "3" /> 
<property name =" VERTICES " > 
<in> 

<v>0 1 2</v> 
<v>0 1 3</v> 
<v>0 2 4</v> 
<v>0 3 4</v> 
<v>l 2 5</v> 
<v>l 3 5</v> 
<v>2 4 5</v> 
<v>3 4 5</v> 
</m> 
</property > 

<property name="RAYS"> 
<m> 

<v>l 0</v> 
<v>0 1 0</v> 
<v>0 l</v> 
<v>0 -l</v> 
<v>0 -1 0</v> 
<v>-l 0</v> 
</m> 
</ property > 

<property name = " N_RAYS " value="6" /> 
<property name=" N.VERTICES " value="8" /> 
< property name=" PARAMETER.POS "> 
<m> 

<v>5 1 -l</v> 
<v>4 2 1 l</v> 
<v>5 2 2 l</v> 
</m> 
</property > 

<property name=" PARAMETER.BOUNDS " > 
<m> 

<v>-12 12</v> 
<v>-12 12</v> 
<v>-12 12</v> 
</m> 
</property > 

<property name =" N.PARAMETER " value="3" /> 
</obj ect > 

Fan 3^4=^62 

<?xml version=" 1 . 0" encoding=" utf -8" ? > 

<object name="f" type =" po lyt ope :: Par amet erFan&lt ; Rat ional &gt ; " version="2 . 9 . 7" 
tm="4b0d42cd" xmlns="http:// www . math. tu-berlin. de/polymake/#3"> 

<property name=" POINT_LIMIT " value="12" /> 

<property name="DIM" value="3" /> 

<property name=" VERTICES " > 
<m> 

<v>0 1 2</v> 
<v>0 1 6</v> 
<v>0 2 5</v> 
<v>0 3 4</v> 
<v>0 3 5</v> 
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<v>0 4 6</v> 
<v>l 2 6</v> 
<v>2 5 6</v> 
<v>3 4 6</v> 
<v>3 5 6</v> 
</m> 
</property > 

<property name="RAYS"> 
<m> 

<v>l 0</v> 
<v>0 1 0</v> 
<v>0 l</v> 
<v>0 -l</v> 
<v>0 1 -l</v> 
<v>0 -1 0</v> 
<v>-l 2 -l</v> 
</m> 
</ property > 

<property name = " N_RAYS " value="7" /> 
<property name = " N_VERTICES " value = "10" /> 
<property name=" PARAMETER.POS " > 
<in> 

<v>5 2 -l</v> 
</m> 
</property > 

<property name=" PARAMETER.BOUNDS " > 

<m> 

<v>-6 5</v> 

</m> 
</property > 

<property name=" N.PARAMETER " value="l" /> 
</obj ect > 
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